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1.  Introduction 


The  traditional  asymptotic  homogenization  method  has  found  wide  acceptance  in  the  composites 
modeling  community  (Bendsoe  and  Kikuchi,  1988;  Fish  et  ah,  1999;  Francfort  and  Murat,  1986; 
Ghosh  et  ah,  1996;  Hollister  et  ah,  1994;  Lene,  1986;  Moulinec  and  Suquet,  1998;  Sanchez- 
Palencia,  1980;  Shkoller  and  Hegemier,  1995),  for  its  ability  to  incorporate  and  estimate 
effective  material  properties  from  a  subscale  representative  unit  cell  (Babuska,  1976).  Central  to 
the  method  is  the  principal  of  convergence,  which  has  been  studied  in  the  mathematical  literature 
for  H-,  G-,  F-,  and  two-scale  convergence.  In  simple  terms,  H-convergence,  for  example,  states 
that  over  an  open  bounded  domain  fl  in  7?^,  a  sequence  of  bounded  and  finite  matrices  is  said 
to  H-converge  to  a  sequence  of  bounded  and  finite  matrices  A  if,  for  all  f  ^  H  (q)  ,  the 
solution  m'"  of  -div{A‘^Vu'^^=  HI{Q.)  converges,  as  the  scaling  parameter  ^ goes  to  zero, 

in  weak-//Q(Q)  to  the  solution  u  of  -div{A'VM)  =  f,ue  and  if  the  sequence  A'^Vu'" 

converges  in  weak-  to  AVu  .  In  general,  the  H-limit  of  A'^ ,  which  is  A,  is  not  explicit 

except  under  certain  assumptions.  The  works  of  Bensoussan  et  al.  (1978)  and  Sanchez-Palencia 
(1980)  are  commonly  cited  to  show  that  the  explicit  H-limit  is  known  for  periodic  systems. 
Loosely  speaking,  for  the  convergence  of  in  the  appropriate  space  of  linear  operators,  or  the 
convergence  of  the  Green’s  operator,  the  literature  refers  to  G-convergence  (Spagnolo,  1976). 

The  convergence  of  [a'^Vu'^ ,  m'")  is  referred  to  as  F-convergence  (Tartar,  1978),  and  a  more 

specific  case  of  H-convergence,  which  takes  into  account  the  periodicity  of  the  solution,  is 
referred  to  as  two-scale  convergence  (Allaire,  1992). 

The  engineering  and  computational  mechanics  literature  have  made  such  extensive  use  of  these 
fundamental  mathematical  developments  that  it  is  difficult  to  provide  a  truly  exhaustive 
bibliography.  We  note  among  many,  however,  the  work  on  structural  optimization  (Bendsoe  and 
Kikuchi,  1988),  biomechanics  (Hollister  et  al.,  1994),  and  multiphysics  solid  mechanics  (Fish  et 
ah,  1999;  Fish  and  Yu,  2001;  Ghosh  et  al.,  1996,  2001;  Lene,  1986;  Moulinec  and  Suquet,  1998; 
Shkoller  and  Hegemier,  1995;  Terada  et  al.,  2000;  Terada  and  Kikuchi,  2001).  Until  now, 
however,  the  engineering  developments  in  homogenization  have  mainly  considered  specialized 
cases  of  two-scale  convergence  where  the  periodicity  is  built  into  the  computational  procedure. 
The  application  of  the  so-called  y-periodicity  generally  precludes  the  ability  to  consider 
nonperiodic  systems.  In  this  regard,  developments  in  the  engineering  and  mathematical  literature 
have  been  on  divergent  paths. 

We  have  found  that  nonperiodicity  can  occur  under  two  scenarios.  The  first  is  when  the  finite 
boundary  has  a  direct  influence  on  the  solution,  or  the  reduction-of-dimension  situation,  which 
can  occur  typically  in  laminated  plates  or  stratified  rods.  The  second  is  when  the  material 
exhibits  finite  microstructure  but  in  a  nonperiodic  sense.  This  can  commonly  be  encountered  in 
the  form  of  localized  behavior,  such  as  that  found  in  damage  mechanics  or  fracture,  or  in 
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smoothly  varying  microscopic  material  dependence  on  macroscopie  coordinates,  such  as  that 
found  in  functionally  graded  eomposites. 

Limited  works  by  the  engineering  community  have  attempted  to  modify  the  homogenization 
method  for  systems  involving  finitely  thiek  plates  (Rostam-Abadi  et  ah,  2000).  Still  others  have 
considered  the  presence  of  strong  anharmonic  local  fields  that  disrupt  the  periodicity  and  require 
increased  resolution  solutions  (Oden  and  Vemaganti,  1999;  Ghosh  et  ah,  2001;  Oden  et  ah, 

2001).  In  these  works,  the  traditional  periodie  homogenization  method  has  been  modified  at  a 
schematic  level  using  the  laminated  plate  theory  (Rostam-Abadi  et  ah,  2000)  or  h  and p  finite 
element  adaptivity  techniques  (Ghosh  et  ah,  2001;  Oden  et  ah,  2001)  so  that  direetions  in 
Cartesian  space  or  regions  within  a  nonperiodically  deforming  body  are  speeifically  aeeounted 
for  by  small  changes  in  the  periodie  homogenization  proeedure.  On  the  other  hand, 
developments  in  the  mathematical  literature  have  been  tied  to  special  cases  where  material 
dependenee  on  spatial  eoordinates  is  strietly  demonstrable  (Briane,  1994;  Shkoller,  1997;  Fabre 
and  Mossino,  1998).  Namely,  they  consider  material  types  where  key  heterogeneous  features  are 
within  a  diffeomorphie  mapping  and  linear  mapping  away  from  a  periodie  arrangement. 

Briane  (1994)  developed  a  nonperiodic  homogenization  method  in  which  a  highly  oscillatory  but 
nonperiodic  matrix  is  compared  to  a  periodie  matrix  B‘'  as  the  mieroseale  feature  size  goes 
(slower)  to  zero  as  s  goes  to  zero.  The  limit  of  is  a  function  of  every  point  in  the  material. 
Recognizing  the  computationally  infeasible  nature  of  this  approach,  Shkoller  (1997)  showed  that 
the  problem  eould  be  simplified  to  a  unit  eell  approaeh,  where  the  mieroseale  feature  size  is  kept 
finite  and  the  error  in  the  approximation  is  of  the  order  of  the  unit  cell  size.  In  this  sense,  the 
asymptotie  nature  of  homogenization  is  not  being  fully  assumed.  Fabre  and  Mossino  (1998) 
eompleted  a  similar  study,  looking  at  H-eonvergent  multiplieable  matriees  that  may  be  used  to 
perform  a  linear  mapping  that  takes  a  periodic  material  into  a  nonperiodic  one  thereby  allowing 
periodic  homogenization  of  a  mapped  nonperiodic  material.  Gustafsson  and  Mossino  (2003) 
developed  a  nonperiodic  homogenization  method  for  diffusion  equations  based  on  the  H- 
convergent  multiplieable  matrices  of  Fabre  and  Mossino  (1998)  to  write  explicit  expressions  of 
homogenization  for  a  plate  and  less  explieit  expressions  for  a  thin  eylinder. 

Developments  in  nonperiodic  homogenization  stem  from  the  general  H-convergence  properties 
of  the  homogenization  theory  eombined  with  mapping  principles  and  dimension  reduetion.  In  its 
most  general  and  primitive  form,  homogenization  makes  no  presumption  of  periodicity,  which 
indicates  that  its  application  to  nonperiodic  systems  can  be  made  forthwith.  Bensoussan  et  al. 
(1978)  have  shown  a  one-dimensional  result  for  homogenization  that  is  implicit  in  its  absence  of 
periodicity  constraints  on  the  microstrueture. 

In  this  report,  we  develop  a  one-dimensional  nonperiodie  computational  homogenization  method 
with  explicit  equations  and  show  that  it  recovers  the  periodic  case.  This  is  done  by  relaxing  the 
periodicity  assumption  of  the  mieroseale  displaeement  field,  integrating  analytically,  and 
applying  a  boundedness  argument  to  the  result.  The  procedure  is  entirely  different  from  that  of 
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Bensoussan  et  al.  (1978)  yet  is  shown  to  yield  the  same  result.  This  enables  the  development  of 
the  new  computational  homogenization  procedure  for  the  one-dimensional  continua  (i.e.,  for  the 
elastic  rod).  Generalization  to  the  three-dimensional  continua  is  expected  to  follow  in  a  later 
paper.  When  compared  to  the  engineering  literature,  the  procedure  described  differs  in  two 
ways:  (1)  the  approximation  space  is  introduced  early  in  the  procedure,  and  (2)  no  periodicity 
condition  is  imposed  on  the  microscale  displacement  field.  This  second  difference  is  expected  to 
make  the  new  procedure  relevant  to  a  wider  class  of  problems  than  the  preceding 
homogenization  procedures.  Using  the  example  of  an  end-displaced  rod  with  linearly  varying 
axial  rigidity  EA,  a  comparison  is  made  between  the  displacement  field  predicted  by  the  new 
homogenization  procedure  and  the  displacement  field  of  the  exact  solution.  A  second  example 
of  a  rod  with  periodically  varying  axial  rigidity  is  used  to  discuss  the  appropriateness  of  the 
microscale  periodicity  condition  imposed  by  preceding  homogenization  procedures. 


2.  On  a  New  Homogenization  Procedure 


2,1  Formulation 

Consider  a  straight  rod  characterized  by  an  axial  rigidity  EA,  whose  deformation  is  described  by 
a  displacement  field  u'^{x‘^)  governed  by  the  equilibrium  equation 


a 

dx^ 


( 

EA 

V 


=  -/, 


(1) 


where  for  the  simplicity  of  discussion,  the  external  forcing  /  will  be  considered  to  be  absent. 
Consider  further  the  rod  to  be  fixed  rigidly  at  its  initiating  end  with  a  prescribed  displacement  at 
its  terminating  end: 

M^(x^=0)  =  0,  (2) 


and 


u‘'(x'^  =  l)  =  u,. 


(3) 


This  choice  of  end  conditions  is  motivated  by  an  eventual  finite  element  implementation  of  the 
described  procedure,  in  which  the  end  conditions  take  the  form  of  equations  2  and  3,  and  a  rigid 
body  translation.  The  axial  rigidity  EA  of  the  described  rod  is  considered  to  be  nonuniform.  It 
can  be  decomposed  into  a  constant  macroscale  axial  rigidity  and  a  nonconstant  microscale  axial 
rigidity 

EA  =  EA^  +  EA\y) ,  (4) 
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where  dependenee  on  the  maeroseale  implies  x  =  x'"  dependenee  and  dependenee  on  the 
mieroseale  implies  y  =x^  j a  dependenee,  ^  being  a  small  parameter.  The  equilibrium  equation 
ean  be  updated  to  refleet  dependenee  on  the  two  seales  by  using  the  derivative  expression 

d  _  d  ^  \  d 
dx‘^  dx  a  dy 


As  in  previous  homogenization  proeedures,  the  displaeement  field  is  expanded  using  a:. 

u‘^{x‘^)  =  u^{x,y)  + au^{x,y)^ — ,  (6) 


where  u^{x,y)  and  u{x,y)  denote  maero- and  mieroseale  displaeement  fields,  respeetively. 
Interehangeably,  u^{x,y)  will  be  termed  the  O(a^)  displaeement,  where  0( )  denotes  the  order 
operator,  and  u\x,y)  will  be  termed  the  O(a^)  displaeement.  Introdueing  equations  5  and  6  into 
equation  1  while  negleeting  quantities  of  O(^)  in  equation  6,  followed  by  a  separation  of  orders 
of  a,  yields  the  following  governing  equations: 


for  a 


dy 


EA 


du 

dy 


=  0; 


(V) 


for  a  ^ 


dEA 


dy 


du^  du 
+ 


A 


dx  dy 


2..0 


+  2EA 


d^u 

dxdy 


dy^ 


■  +  EA 


=  0; 


(8) 


for  a 


and  for  a\ 


dEA  du 
dy  dx 


+  EA 


a2  0 

d  u 


+  2EA 


a2  1 

d  u 
dxdy 


=  0; 


EA 


d^u' 


=  0. 


(9) 


(10) 


Without  loss  of  generality,  the  boundary  inhomogeneity  u,  is  not  expanded  and  the  boundary 
eonditions  ean  be  restated  as 

m“(0,0)  =  0, 

Myo.o;  =  o, 

u^QJj a)  =  u^i , 


and 
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(11) 

(12) 

(13) 


u{l,lls)  =  0 . 


(14) 


Integrating  equation  7  twice  with  respect  to  y  yields 

M“(x,y)  =  Ci(x)«(y)  +  C2(x),  (15) 

where  a(y)  =\—dy.  Assuming  C,  (x)  is  not  proportional  to  s,  so  as  to  preserve  the  relative 

J  EA 

ordering  of  all  terms  in  (equation  6)  for  ally,  the  presence  of  a{y)  in  equation  15  results  in 
M°(x,y)  being  unbounded  as  y  ^  oo  ,  or  equivalently  then,  as  ^  0  .  For  example,  retaining 
only  the  macroscale  term  in  the  axial  rigidity  for  ease  of  explanation,  at  the  terminal  end 

a(y  =  lls)  evaluates  to  which  is  unbounded  as  ^  0  .  We  require  that  the  solution 

procedure  be  bounded,  or  appealing  to  a  physical  argument,  be  applicable  to  the  case  of  a 
homogenous  rod  where  ^is  arbitrarily  small,  in  which  case,  the  need  for  a  bounded  M°(x,y) 
requires  that  Cj(x)  be  zero. 

To  further  facilitate  a  solution,  an  approximating  space  is  introduced  for  C2(x) .  Specifically,  we 
search  for  an  approximation  to  M°(x,y)  such  that  C2(x)  e  //'(x)  (i.e.,  the  Sobolev  space  of 
degree  1)  (Hughes,  1987).  For  simplicity  herein,  a  single  element  that  spans  the  rod’s  length  / 
will  be  considered,  but,  more  generally,  the  use  of  the  77 '  Sobolev  approximation  space  admits 
usage  of  multiple  elements  while  insuring  displacement  continuity  across  the  element 
boundaries.  In  accordance  with  this  discussion,  a  linear  polynomial  approximation  for  C2(x)  is 
introduced: 


C2(x)  =  a2X  +  ^2, 


(16) 


where  a2  and  62  are  unknown  constants.  Applying  the  boundary  conditions  (equations  1 1  and 
13)  yields 

u\x,y)  =  u\x)  =  ^x.  (17) 


With  the  determination  of  the  O(s^)  displacement  completed,  the  analysis  of  the  O(s^) 
displacement  follows.  Substituting  equation  17  into  equations  8,  9,  and  10  yields 


dEA 


7..0 


dy 


Ml  du 


I  dy 


^  +  EA^-^  =  Q, 


dy^ 


dEA  du'  ^  ^  .  d^u' 
dy  dx  dxdy 


0, 


(18) 

(19) 


and 


a2  1 

£.4^  =  0, 

5x' 


(20) 
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Equation  1 8  can  be  rewritten  as 


A 

dy 


( 

EA 

V 


du  ^ 


d  ( EA 
—  — u 

/ 


0 

i 


\ 

y 


which  can  be  integrated  once  and  rearranged  to  yield 

du\  EA  0 

EA  —  +  —  U,  =C^ix). 
dy  I 


(21) 


(22) 


Note  that  equation  22  states  that  a  function  of  x  andy  (say  g(x,y)  =  EA—+—u°)  is  equal  to  a 

dy  I 

function  of  x  only  (say  f(x)  =  C-^{x) )  for  all  x,y — this  can  imply  that  either  g(x,  yj  is  a  funetion  of 

X  only  or  both  gfx,  y)  and  f(x)  are  constant.  The  former  does  not  lead  to  an  immediate  solution,* 
while  the  latter  leads  to  the  solution 


C3  (x)  =  C3 ,  u '  (X,  y )  =  C3  «(y )  -  ^  u °  +  Q  (x)  ,  (23) 

where  C3  is  a  pure  constant  to  be  determined  by  the  boundary  conditions.  The  substitution  of 
equation  23  into  equation  19  and  the  enforcement  of  the  resulting  equation  for  ally  determine 
that  C4(x)  =  C4  (i.e.,  that  C4  is  also  a  pure  constant  to  be  determined  by  the  boundary 
conditions).'!'  Note  that  the  resulting  solution  for  the  microscale  displacement  u  is  independent 
of  X  and  satisfies  equation  20  identically.  Following  application  of  equations  12  and  14,  the 
microscale  displacement  is  given  by 

sa{lls)  I 

Due  to  the  presence  of  the  multiplier  s  in  equation  6,  s  m'  (y)  is  bounded  for  all  y  as  ^  0  .  Note 
that  the  microscale  displacement  is  not  necessarily  periodic  iny  even  if  the  axial  rigidity  EA  is. 
Traditional  homogenization  procedures  summarily  impose  periodicity  on  the  microscale 
displacement.  As  a  final  note,  for  the  ease  of  a  homogenous  rod  (i.e.,  EA'^{y  )  =  0),  equation  24 
evaluates  to  zero,  the  macroscale  displacement  is  given  by  equation  17,  and  the  exact  solution  is 
recovered. 


2,2  Extension  of  Procedure  to  Include  Higher-Order  Terms  in  u^(x^ ) 

The  solutions  presented  in  seetion  2.1  were  found  by  truncating  the  asymptotie  expansion 
(equation  6)  up  to  O(^).  If  instead  truncation  is  not  performed,  a  new  term  involving  m^(x,  y) 


Simplification  of  the  equations  does  not  occur  for  this  choice,  but  the  possibility  of  a  second  solution  arising  from  this  choice 
may  exist. 

^For  the  case  of  a  nonzero  external  \oa.Af(x),  C.  is  replaced  by  ^  ^  _ i_  (y  . 

dEA/dy  “ 
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appears  in  equations  9  and  19.  Similarly,  new  terms  involving  u^{x,y)  and  u^{x,y)  appear  in 
equations  10  and  20.  Note  that  equation  10,  and  thus  equation  20,  were  not  used  in  developing 
the  solution  proeedure  but  instead  were  eheeked  for  satisfaetion.  In  the  solution  presented  for 
this  seetion,  attention  is  therefore  focused  only  on  solving  the  new  equation  19,  namely 


BE  A  5m' 
dy  dx 


+  2EA - +  — 

dxdy  dy 


EA 


du 


2  A 


dy 


=  0, 


(25) 


with  the  still  valid  equations  22  and  23.  In  this  manner,  macro-  and  microscale  displacements 
will  again  be  developed  without  the  need  to  first  truncate  the  displacement  asymptotic  series 
given  by  equation  6. 


Substituting  equation  23  into  equation  25  and  redistributing  yields 


C  = 

'-'4 


dy 


V  J 


dEA 

dy 


(26) 


where  a  prime  denotes  differentiation  with  respect  to  x.  Note  that  equation  26  states  that  a  function 
ofx  (say  f(x)  =  C4)  is  equal  to  a  second  function  of  both  x  andy  (say  g(x,y)  =  -AfEA^V^)- 

ey(  dy  )l  dy 

As  in  the  previous  section,  a  solution  in  which  both  f(x)  and  g(x,y)  are  the  same  constant  is  sought. 
This  leads  to  the  expression  C4  =  C^x  +  Cg  and  an  updated  expression  for  M'(x,y) , 


u\x,y)  =  C^a{y)-jU^  +C^x  +  Cf,, 


(27) 


with  three  undetermined  true  constants.  Application  of  the  first  boundary  condition 
(equation  12)  yields  Ce  =  0,  while  the  second  boundary  condition  (equation  13)  yields 


C3  a 


rn 

UJ 


-uf+CJ  =  0. 
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(28) 


Separating  equation  28  into  0(s  ')  and  0(^')  terms  yields  two  equations  for  the  remaining 
constants: 


for  s  ', 


and  for  £■", 


C3  a\ 


rn 

UJ 


1  0  ^ 

—  M,  =  U  ; 
a 


C5/  =  0 , 


(29) 


(30) 
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which,  when  satisfied  again,  yields  equation  24. 

2,3  Application  to  an  Example  Rod — Linearly  Varying  Axial  Rigidity  EA 

The  new  homogenization  proeedure  deseribed  in  seetion  2.1  is  applied  to  a  rod  of  unit  length 
with  a  linearly  inereasing  axial  rigidity  given  by 

£^  =  3  +  0.01y,  (31) 

and  with  an  end  displaeement  given  by  =  0.1 .  The  microsealey  is  chosen  to  be  one  one- 
hundredth  of  the  maeroscale  x  such  that  £■  =  0.01 .  This  information,  evaluated  with  the 
definition  of  a(y)  ,  the  truncated  expansion  (equation  6),  and  the  ordered  displacements 
(equations  17  and  24),  yields  a  predieted  solution  for  the  displaeement  field 

=  -0.382  +  0.348 ln(3  +  ,  (32) 


where  all  ealculations  have  been  earned  out  to  three  signifieant  figures.  For  comparison,  an 
exact  solution  exists  for  the  linearly  varying  rod.  Colleeting  the  relevant  information,  the 
boundary  value  problem  governing  the  rod’s  exact  displacement  field  is  as  follows; 


dx 


EA  — 
dx 


=  0,  EA  =  3  +  0.01y,  u(0)  =  0,  m(1)  =  0.1, 


(33) 


which  has  the  solution 


1  ln(3)  1  ln(3  +  x) 

To  ln(3)-21n(2)  To  ln(3)  -  2  ln(2)  ' 


(34) 


Evaluating  equation  34  to  three  significant  figures  yields  the  exact  displacement 

m(x)  = -0.382  + 0.348  ln(3  +  x),  (35) 


whieh  matehes  the  perturbation  solution  (equation  32).  The  same  mateh  oecurs  when  the 
ealeulations  are  earned  out  to  any  number  of  signifieant  figures,  indieating  that  for  the  linearly 
varying  rod,  the  newly  described  homogenization  proeedure  reeovers  the  exact  solution. 

2,4  Application  to  an  Example  Rod — Periodically  Varying  Axial  Rigidity  EA 

As  a  second  example,  the  new  homogenization  procedure  is  applied  to  a  rod  of  unit  length  with  a 
periodically  varying  axial  rigidity  given  by 

E^  =  3  +  sin(^),  (36) 

and  with  an  end  displacement  given  by  =  0.1 .  The  microscale y  is  chosen  to  be  one  one- 
hundredth  of  the  maeroseale  x  sueh  that  ^  =  0.0 1  (note  that  this  gives  a  y-period  of  1 0  for  the 
axial  rigidity).  This  information,  evaluated  with  the  definition  of  «(y)  ,  the  truncated  expansion 
(equation  6),  and  the  ordered  displacements  (equations  17  and  24),  yields  a  predieted  solution  for 
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Figure  1.  (a)  The  function  a{y),  (b)  the  microscale  displacement  field  su^iy),  and  (c)  the  predicted  displacement 
field  uXx^^  for  the  example  rod  with  periodically  varying  axial  rigidity  (i)  EA  =  3  +  sin(;ry/5^, 

(ii)  EA^3  +  2.9sin(;r>’/5),  and  (iii)  _£4  =  3  +  2.9signum(sm(n:y/5y). 

the  displacement  field.  The  top  portion  of  figure  1  shows  the  results  for  the  function  a{y)  ,  the 
predicted  microscale  displacement  field  su^iy),  and  the  predicted  displacement  field  uXx‘X  for 
this  example. 

ft  is  interesting  to  note  from  this  example  that  the  microscale  displacement  is  indeed  jp-periodic, 
as  imposed  a  priori  in  the  traditional  homogenization  procedure.  The  effects  of  further 
increasing  the  variability  of  the  axial  rigidity  are  shown  in  the  middle  portion  of  the  figure. 
Finally,  the  results  for  a  rod  with  nonsinusoidal,  but  still  periodic,  axial  rigidity  are  given  in  the 
bottom  portion  of  the  figure.  This  last  example  gives  further  evidence  supporting  the  imposition 
of  microscale  displacement  periodicity  in  the  standard  homogenization  method  (note  that, 
indeed,  the  microscale  displacement  is  periodic). 
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3.  Conclusions 


Using  the  well-established  homogenization  theory,  a  modification  has  been  proposed  for 
engineering  applications  in  which  the  microscale  displacement  field  is  nonperiodic.  The  result  is 
a  novel  nonperiodic  homogenization  procedure  for  functionally  graded  and  reduced  dimension 
problems  where  y-periodicity  breaks  down  because  of  material  or  dimensional  considerations.  A 
specialized  one-dimensional  formulation  was  presented,  verified  against  an  exact  solution,  and 
subsequently  demonstrated.  A  fully  three-dimensional  formulation  is  currently  underway  and 
will  be  presented  in  a  sequel  paper. 
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